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"q^" ■ Abstract 
ON _ 

The analysis of the statistical and dynamical fluctuations in nucleus-nucleus 
collisions on an event-by-event basis strongly relies on a comparison with spe- 
cially constructed artificial events where statistical fluctuations and kinemat- 
ON ' ical correlations are under control. In this paper, we present a novel, analyt- 



(N 



ical method of constructing reference events based on independent emission, 
. modified by the energy /momentum constraint, which can lead to a better 

understanding of the nature of the observed final-state fluctuations. This ap- 
Q proach can be easily used in the analysis of other topics in the heavy ion field 

' (e.g. flow, HBT etc.) allowing more precise measurements. 

ON 

An event-by-event analysis (e-b-e), so successfully used from the very beginning of high 
^ | energy physics (bubble and streamer chamber experiments), was recently proposed [0J and 
applied B to ultra-relativistic heavy nucleus-nucleus collisions at the CERN SPS energy. 
In hadronic physics, the e-b-e results provided information on the properties of individual 
interactions and their variations from event to event. In the heavy ion field, the search 
for "unusual" events, i.e. events having a particularly high variation of some variable from 
its average value, is especially important due to the expectation of non-trivial dynamical 
fluctuations caused by the formation of Quark-Gluon Plasma (QGP) bubbles and/or other 
exotic phenomena, such as a disoriented chiral condensate (DCC) ||, jet quenching [|J, 
color fluctuations in the early stages of the collision || and others. Recent data from the 
CERN SPS indicate that the energy density reached in nucleus-nucleus collisions has already 
exceeded the estimated critical value for QGP formation in the "average" event, thus the 
e-b-e approach may allow one to separate events in which plasma was created from those in 
which it was not. Futhermore, it may allow one to determine to what degree the transient 
QGP phase and the dynamics of the phase transition to hadronic matter affect fluctuations 
observed in the hadronic state at freeze-out. 

The obvious, necessary experimental condition for e-b-e analysis is a large phase-space 
acceptance. However, even having a large fraction of the charged particles under control, 
one is confronted with the uncertainties related to the neutral, usually undetected, particles 
produced in the collision. Therefore, it is essential to determine to what degree the fluc- 
tuations present in the experimental data (from the measurable part of the phase space) 
represent those of the entire event (relevant for comparison with the theoretical predictions). 
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In this paper we outline a novel, analytical method of deriving complete information 
on the entire reference (so called "mixed") event, based exclusively on the available experi- 
mental information - the fraction of final-state particles observed in the detectors - and the 
central limit theorem. Our approach, as described below, substantially refines the traditional 
one, where mixed events were composed by drawing particles randomly from the huge pool 
created by combining large amount of the data from the same trigger. Most importantly, 
we impose constraints from the conservation laws with associated kinematical correlations; 
this to the best of our knowledge, has never been done before. Some aspects of this concept 
have been already applied in our earlier work (flow analysis of the Bevalac streamer chamber 
data) §. 

In order to establish whether the observed fluctuations are partly dynamical in nature, 
we need to disentangle statistical effects i.e. effects due to the finite number of particles 
in the final state of the collision. In the following we will concentrate on the transverse 
momentum distribution to demonstrate the method qualitatively; it can also be applied to 
other observables. 

First, we define a scale for measuring fluctuations and subtracting the trivial, statistical 
effects from the overall event-to-event variations. We start by comparing the width of the 
experimental spectrum with that of the specially constructed, mixed events of the same 
multiplicity for which we assumed independent particle emission, modified by the momen- 
tum/energy (pj;,pj,,p z ,E) conservation laws[]. In mixed events, each particle is sampled from 
a different heavy ion collision belonging to the same data set (data taken with the same 
t rigger )0. 

N 

i=i 

where 

pmix _ density Q f a }} qsj) particles of the mixed event 
p 1 , p 2 , ... - particle momenta 

p\ (ft) - density of a single particle having momentum p, in the i-th event 
5(P — Eft) - imposes momentum/energy conservation 

P - total momentum/energy of the initial state = Pbeam in a fixed-target experiment 

Eft- total momentum of all particles in the final state. 
While dynamical fluctuations were totally eliminated from the mixed events during construc- 
tion by the independent emission mechanism, we carefully preserved their exact multiplicity: 
i.e., for each data event of multiplicity N, we took one particle from N different data events 
with the same trigger to compose an equivalent mixed event of N particles. 



lr rhe momentum/energy conservation explicitly addresses correlations particularly related to the 
observable pt chosen for our example. To investigate different observables, other constraints need 
to be included: e.g. the K/ir ratio requires the imposition of strangeness conservation. 

2 If N is large - mixed and real (data) events have the same single-particle spectra; if N is small - 
one needs to apply a correction factor of order (N-l)/N. 
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The width of the event distribution constructed thusly provides a benchmark for our 
analysis. It represents the scale of the statistical fluctuations and a measure of the probability 
for events to appear in the tails of the spectrum (rare events), where the fluctuations are 
largest. 

The most important advantages of this strategy include: 

• model independence 

• the same experimental systematics (efficiencies, acceptances, resolutions, etc.) in data 
and mixed events 

• consistency with the conservation laws (via the 5 function), not only excluding "non- 
physical" events from the mixed-events sample, but also introducing into mixed events 
the kinematical correlations naturally present in the data. 

We use the following notation: 

each event consists of N particles 

1, 2, 3, m denote the charged particles observed in the experiment 
m+1, m+2, N denote the non-observed particles (neutral, outside the acceptance, 
etc.) 

We express the density of observed particles as the density of all particles integrated over 
the unobserved part of phase space using the central limit theorem (for simplicity we drop 
the superscript "mix" from the left-hand side): 

Pm(Pi,P2,-Pj = J P N (PnP2T-P m iP m+ n --Pm) d P m+ i-- d P N = 

Next, we replace the density of all particles p N by the superposition of independent single 
particle densities (assuming independent emission) and separate observed and unobserved 
particles. We also write down momentum/energy conservation explicitly. 

/m N m JV 

T[p k M II p k M) s ( p -J2^~ Hpi) d p m+ i-dp N = 



i — 1 m + 1 



The product of the densities of all unobserved particles J f p\ (p k ) is substituted by one ran- 

k — m+ 1 

N 

dom variable p(E pu) representing the total missing momentum in the collision. Note, that 

m+1 

/JV „ N 

H P k i.P k )dp--- = / p(E^ p fe )d(Ep fe ). The new variable p(^Pk) reflects much better 

k — m + 1 m+1 

the measuring capabilities of the experiment. 

Multiplicity at CERN SPS energies is high, ranging up to a few thousand particles in 
central Pb+Pb collisions; therefore, k is large enough to allow us to use the central limit 
theorem to approximate the density of unobserved particles by a normal distribution Af: 

JV JV 

p( E Pk) ~A/"( ]T Pk ) 
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where 

M{y) = p . e -^ A ^- {ymyk - {y)k) . 
A ik is the inverse of the covariance matrix 

A = C' 1 

C = cov(y it y k ) = (y l - (yY)(y k - (y) k ) 
and i, k = 1, 2, 3, 4 (e.g., p x ,p y ,p z , E). 

We assumed that all errors are the same, therefore, all weights are the same. 

m n N m N N 

= n rip*) I W E pMp - Ep* - E Pk) d( E p*) = 

k — 1 k — TTi+1 k — 1 fc— m + 1 fc— m+1 

after integrating with the 5 function, we obtain 

m m 

= rip fc (pj^-E^)- 

Note, that by using the central limit theorem and by integrating with the 5 function we 
were able to eliminate all variables related to experimentally inaccessible particles. The final 
results: 

m m 

Pm(Pl,P2, -Pm) = II P'iPk) N( P ~ T,Pk) 
fc=l fc=l 

depend only on the single-particle quantities measured in the experiment and on the weight- 
ing factor, also totally calculable from experimental data. Thus, our task narrows down to 

m 

computing the weighting factor W = M{P — ^pt). 

fc=i 

Before we discuss the practical aspects of computing W, let us make two digressions on 
statistics: 

• W has a number of very convenient features: 

-We (0,1) 

— it eliminates the events where momentum/energy are poorly conserved 

— it treats properly correlations resulting from conservation laws 

— it is totally known from the experiment 

• The covariance matrices for the sum over the observed particles and for the sum over 
the unobserved particles are equal. 
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Having sketched a general outline of the new method, let us point out the necessary steps 
to calculate weights (W) and to construct proper reference events. 
The entire procedure is factorized to the five steps: 

m 

step 1: For all observable particles in each event calculate W = J^Pfc where each p^. is a 

n=l 

momentum-energy four vector. 

step 2: Calculate the averages (W l ) and elements of the covariant matrix for the entire 
set of event s0: 

c ik = ((W* - (W*))(W k - (W k ))) 

where W is a four vector and i,k denote its components. 

step 3: Calculate the inverse of the covariance matrix: A = C 1 (A-C = I). 

step 4: Find W (= A/"(W) for each event 

M = const • e -¥w-(w)rA(w- { w)) ( W = four yector) 

to be used in: 

in 

p m (Pi,p>, -pj = n p k (p k ) ■ wo. 
*,=! 

step 5: The last step we call "simulations of measurements": Let us assume that our cal- 
culations show that some particular "mixed" event has the probability A/"(W) of 0.3. We 
need to convert this number to 1 (= event entering our reference data sample) or (= event 
rejected) to have a uniform treatment with the experiment (all events collected on the DST 
have probability = 1; all missing, of course, have probability = 0). This conversion is done 
in the following way: we draw a random number (a) between and 1, and compare our 
calculated weight W=A/"(W) against it. 

If a is smaller than the W of the particular event, this event is accepted with a new 
probability = 1; however if a is bigger than W - the event does not enter our sample of 
reference data. So, in our example, the randomly selected a has to be smaller than 0.3 in 
order for the event to be accepted into the reference data sample. 

The method described above is presently being tested with both Monte Carlo and exper- 
imental data. The quantitative understanding and evaluation of the results will take some 
time. In particular, in order to precisely reproduce the single particle distribution we have 
to introduce an additional weighting factor which compensates for the exponent resulting 
from the integration with the 5 function. However, while still in the process of testing, we 
would like to communicate and make it available to the community, due to its wide range of 
applications for topics other than fluctuations analysis. We expect to present results from 
our simulations and data analysis soon. 
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3 Very schematic. Full procedure is more complicated. 
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